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Q ■ Abstract 

O ' We investigate the integrabihty of the discrete non-hnear equation governing 

^ ' the dependence on geodesic distance of planar graphs with inner vertices of 

. even valences. This equation follows from a bijection between graphs and blos- 

som trees and is expressed in terms of generating functions for random walks. 
We construct explicitly an infinite set of conserved quantities for this equation, 
also involving suitable combinations of random walk generating functions. The 
^ ■ proof of their conservation, i.e. their eventual independence on the geodesic 

■<s:J- . distance, relies on the connection between random walks and heaps of dimers. 

i The values of the conserved quantities are identified with generating functions 

Q I for graphs with fixed numbers of external legs. Alternative equivalent choices 

■ for the set of conserved quantities are also discussed and some applications are 

■ presented. 
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1. Introduction 

Graph combinatorics seems to provide a remarkable source for (discrete) integrable 
systems. So far, the underlying integrability seemed to be intimately related to the exis- 
tence of a matrix model formulation for the counting of these graphs [1,2]. In the context 
of matrix model solutions to 2D Quantum Gravity [3], the generating functions for possibly 
decorated graphs of arbitrary genus may indeed be interpreted as tau-functions of various 
integrable hierarchies [4]. 

More recently, an alternative purely combinatorial and bijective approach to the enu- 
meration of planar graphs was developed, based on the transformation of graphs into 
decorated trees [5-7]. In this context, another remarkable, apparently unrelated integrable 
structure was observed, now involving the geodesic distance on the graphs [8] . 

More precisely, in Ref . [8] , special attention was paid to the case of planar graphs with 
vertices of even valence only, and to their generating function with two distinguished 
points at geodesic distance at most n, with n a non-negative integer. The latter was shown 
to obey a self-consistent master equation turning into a recursion relation on n for graphs 
with bounded valences. This non-linear recursion relation turned out to be exactly solvable 
in terms of discrete soliton-like tau-functions. The precise solution involved integration 
constants which may be rephrased into conserved quantities, thus displaying explicitly 
the integrability of the master equation. The knowledge of these conserved quantities 
provides a powerful tool for generating explicit expressions for various generating functions 
of interest in the graph combinatorial framework. In particular, exploiting the conserved 
quantities allows to bypass the quite tedious use of the exact solutions of Ref. [8] . 

The aim of this paper is to construct ah initio a set of conserved quantities for the 
master equation in a purely combinatorial setting, and to interpret them in the language 
of graphs. As will be recalled below, the master equation has a compact expression in 
terms of discrete random walks. Our construction relies therefore on properties of parti- 
tion functions for such random walks, as well as for "hard dimers" on a line, both with 
suitable weights involving the i?n's. Both objects are known [9,11] to be related via some 
boson/fermion type inversion relation. The origin of this relation is best understood by 
relating random walks to so-called "heaps of dimers" , a boson-like counterpart of the mu- 
tually excluding, thus fermion-like, hard dimers. 

The paper is organized as follows. In Section 2, we first recall the existing bijection 
between on the one hand planar graphs with two external legs and inner vertices of even 
valence, and on the other hand planted trees with two kinds of leaves, the so-called blossom 
trees (Sect.2.1). This bijection is used to re-derive the master equation for in Sect. 2. 2. 
Other possible applications of the master equation are discussed in Sect. 2. 3, as well as its 
integrable nature in Sect. 2. 4. Section 3 is devoted to the study of generating functions 
for random walks, heaps of dimers, and one-dimensional hard dimer configurations and to 
their relationships. In Section 3.1, we first present a useful fundamental relation for the 
generating function of random walks, which is then identified with that of a particular class 
of heaps of dimers called pyramids. Details of this identification are gathered in Appendix 
A. This identification allows to derive a crucial inversion relation involving the partition 
function for hard dimers on a segment. In Section 4, we first give explicit formulas for 
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the conserved quantities T2i{n), i = 1,2,... in terms of partition functions for random 
walks (Sect. 4.1). The actual conservation of these quantities, i.e. their independence of 
n whenever the master equation is satisfied, is proved in Sect. 4. 2 by extensive use of the 
fundamental and inversion formulas of Section 3. Section 5 concerns the graph interpre- 
tation of the conserved quantities T2i{n) in terms of 2z-point functions, i.e. generating 
functions for graphs with 2i external legs (Sect. 5.1). In the case of graphs with bounded 
valences, say up to 2m, the infinite set of conserved quantities is reducible to the finite 
set of the m — 1 first ones. The corresponding (linear) relations between their constant 
values are derived in Sect. 5. 2. Finally, Sect. 5. 3 presents an alternative bijective proof for 
the conservation of r2 (n) , interpreted as the generating function for graphs with the two 
legs in the same face. Possible extensions of this bijective proof to higher order conserved 
quantities are also discussed, in view of particularly suggestive identities equivalent to the 
conservation of for z = 2, 3, 4. In Section 6, we first present two equivalent alternative 
sets of conserved quantities, one respecting the "time-reversal" (n — > — n) symmetry of the 
master equation (Sect. 6.1) and the other compacted, i.e. involving the least possible num- 
ber of RnS (Sect. 6. 2). Proofs of their conservation rely on generalized inversion formulas, 
all displayed in Appendix B. To illustrate the interest of using conserved quantities, we 
give in Sect. 6. 3 an example of application to the statistics of neighbors of the external face 
in pure tetravalent and pure hexavalent graphs. We gather a few concluding remarks in 
Section 7. 

2. Generating functions for planar graphs: master equation 

2.1. From planar graphs to blossom trees 

Let us first recall the bijection between planar graphs and blossom trees. More pre- 
cisely, by planar graphs, we mean here graphs with planar topology, with inner vertices of 
even valences only and with two univalent endpoints (see Fig.l-(a) for an example). We 
shall refer to the endpoints as legs and to the graphs as two-leg diagrams. The two legs 
are distinguished into an incoming and an outcoming leg and need not be adjacent to the 
same face. We call geodesic distance between the two legs the minimal number of edges 
to be crossed in a path joining both legs in the plane. In the planar representation, we 
choose to have the incoming leg adjacent to the external face. 

On the other hand, by blossom trees, we mean planted plane trees satisfying 

(Bl) all the inner vertices of the tree have even valences; 

(B2) the endpoints of the tree are of two types: buds and leaves; 

(B3) each 2A;-valent inner vertex is adjacent to exactly k — 1 buds; 
(see Fig.l-(e) for an example). In this characterization, it is assumed that the tree has at 
least one inner vertex but, by convention, the tree reduced to a single leaf attached to the 
root is also considered as a blossom tree. Note that from property (B3), we deduce that 
the total number of leaves in a blossom tree is equal to the total number of buds plus one. 

The bijection may be obtained by cutting each two-leg diagram into a blossom tree as 
follows. We first visit successively each edge bordering the external face in counterclockwise 
direction around the diagram (see Fig.l-(b)), starting from the incoming leg. Any visited 
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(a) (b) (c) 




(d) (e) 

Fig. 1: A two leg-diagram with distinguished incoming and outcoming legs 
(a). The incoming leg is adjacent to the external face. All inner vertices 
have even valence. Here the geodesic distance between the legs is 1. This 
diagram is transformed into a blossom tree (e) by cutting edges into bud-leaf 
pairs as explained in the text. Buds (resp. leaves) are represented by black 
(resp. white) arrows. Here the cutting procedure requires performing two 
turns around the graph: (b)^(c) and (c)— ^(d). The outcoming leg serves as 
root for the blossom tree while the incoming one is replaced by a leaf. In the 
blossom tree, each inner vertex of valence 2k is adjacent to exactly k — 1 buds. 

edge is cut into a bud-leaf pair if and only if this operation does not disconnect the diagram. 
After one turn, the net result has been to merge a number of internal faces with the external 
one. We repeat the process on this newly obtained diagram (see Fig.l-(c)) until all faces 
are merged, resulting into a tree (see Fig.l-(d)). We finally replace the incoming leg by 
a leaf and plant the tree at the outcoming one (see Fig.l-(e)). As shown in Ref.[5], the 
resulting tree is a blossom tree and the above construction establishes a bijection between, 
on the one hand two-leg diagrams and, on the other hand blossom trees. Note that the 
inner vertices of the obtained tree are clearly in one-to-one correspondence with those of 
the corresponding diagram and that their valence is preserved. This allows to keep track 
of vertex valences in both pictures. 
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(c) (d) 

Fig. 2: The closing procedure from a blossom tree (a) to a two-leg diagram 
(d). Each pair of a bud immediately followed counterclockwise by a leaf is 
glued into an edge (b). This is iterated (c) until all leaves but one are glued. 
This remaining leaf is replaced by an incoming leg while the root becomes an 
outcoming one. The closing procedure here encircles the root once, hence the 
geodesic distance between the legs is 1. 

The inverse of the cutting procedure consists in reading the sequence of buds and 
leaves in counterclockwise order around the tree (see Fig.2-(a)) and connecting each pair 
of a bud immediately followed by a leaf into an edge (see Fig.2-(b)) The process is repeated 
until only one leaf is left (see Fig.2-(c)), which is chosen as the incoming leg while the root 
becomes the outcoming one (see Fig.2-(d)). Note that this root is usually encircled by a 
number of edges separating it from the outer face. This number is nothing but the geodesic 
distance between the legs. A simple way to measure it is to forbid to encircle the root in 
the closing process. The geodesic distance is then given by the number of excess buds, 
i.e. those buds which remain unmatched. Alternatively, we may define the contour walk 
of a blossom tree by reading the sequence of buds and leaves clockwise around the tree 
starting from the root and performing a down (resp. up) step for each encountered bud 
(resp. leaf) (see Fig. 3). Starting from height and making steps of ±1, the contour walk 
ends at height 1 and the geodesic distance between the legs of the corresponding diagram 
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height > -1 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ depth < n 

Fig. 3: The contour and surrounding wallcs of tlie blossom tree of Fig.l(e), 
or equivalently of Fig. 2(a). The contour walk, represented in solid line, is 
here to be read from right to left and records the succession of buds and 
leaves clockwise around the tree from the root. Starting from height 0, it 
makes an ascending step for each encountered leaf and a descending one for 
each bud, hence ends at height 1. Its depth is here 1, equal to the absolute 
value of the minimum visited height. The surrounding walk, represented in 
dotted line, is to be read from left to right and records to sequence of buds 
and blossom subtrees counterclockwise around the root vertex from the root. 
To check whether the tree at hand contributes to we let the surrounding 
walk start from height n and make an ascending step for each encountered 
bud and a descending one for each blossom subtree. Each descending step 
may be viewed as the net contribution of the subtree within the contour walk. 
Here each blob corresponds to such a subtree. The depth of the contour walk 
will remain less or equal to n if the depth of each (vertically shifted) contour 
walk in a blob starting at height i remains less or equal to i. 



is given by the depth of the walk, defined as the absolute value of the minimum height 
reached by the walk. 



2.2. Master equation 

We may now define the generating function Rn for two-leg diagrams with weight gk 
per 2/c-valent vertex and with legs at a geodesic distance less or equal to n. From the 
above bijection, Rn is also the generating function for blossom trees with weight gk per 
2/c-valent inner vertex and whose contour walk has depth less or equal to n or equivalently 
with at most n excess buds in the closing process that forbids encircling the root. 

It is now easy to write down a master equation for R^ by decomposing each blossom 
tree according to the environment of the inner vertex attached to the root (root vertex). 
This environment may be decomposed into a sequence of buds and descendent subtrees. 
All these subtrees are themselves blossom trees either made of a single leaf or satisfying 
the above characterizations (B1-B3). Note that the contour walk of the original tree is the 
concatenation of down steps for the buds and of the (vertically shifted) contour walks of 
the blossom subtrees. Hence, the constraint of depth less than n translates into constraints 
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on the depth reached by the contour walk of each blossom subtree. This turns into a closed 
relation between the -R^'s- 

More precisely, the environment of the root vertex is best described by yet another 
walk, the surrounding walk, now obtained by reading the sequence of buds and subtrees 
in the opposite, counterclockwise direction around the root vertex. Starting from height 
n at the root, we make a step +1 for each encountered bud and —1 for each encountered 
blossom subtree (see Fig. 3). Each —1 step may be seen as the net contribution of the 
blossom subtree at hand within the total (reversed) contour walk. To ensure that the 
depth of the entire contour walk remains less or equal to n, we must require that each 
subtree encountered at height i in the surrounding walk has a contour walk with depth 
less or equal to i. This leads us to attach a weight Ri to each —1 step from height i to 
height (i — 1). 

The description of the root vertex environment by its surrounding walk and the at- 
tached weights are efficiently encoded into a "transfer matrix" Q acting on a formal or- 
thonormal basis \i) {i E 2) as 

Q\i) = \i+l)+R.\i-l) (2.1) 

The generating function Za,b{k) for walks of k steps going from, say height a to height b, 
and with weight Ri for each descent z — > (^ — 1) is easily expressed in terms of Q as 

Za,b{k) = {blQ'^la) (2.2) 

Note that Za,b{k) is non zero only for {b — a) = k mod 2. From the above analysis, the 
quantity gkZn,n-i{'^k — 1) is the generating function for blossom trees with a 2A;-valent 
root vertex and with contour walk of depth less or equal to n. To incorporate arbitrary 
(even) valences, we introduce the quantities 

V'{Q)^Y.9kQ^''~^ (2.3) 

and 

<6 = Wma) = Y,9kZa,b{2k - 1) (2.4) 

fc>i 

(which is non zero only for (6 — a) odd). The quantity V^^ may be viewed as a "grand- 
canonical" generating function for walks from a to 6 and with arbitrary odd length. 
With these notations, we may finally write the master equation 



Rn-l + 



(2.5) 



where the first contribution (1) corresponds to having a single leaf connected to the root. 
This relation is valid for all n > with the convention that i?j = for alH < 0. This 
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Fig. 4: A schematic representation of the truncated master equation (2.6). 
The generating function Rn is decomposed according to the environment of 
the root vertex, with 0, 1, 2 buds and 1, 2, 3 blossom subtrees according 
to its valence 2, 4, 6 respectively. To ensure that the depth of the contour 
walk does not exceed n, we have to impose bounds on the maximal number 
of excess buds in each subtree. We have represented (here for n = 3) these 
excess buds in maximal number for each subtree. The index of R for each 
subtree is nothing but this maximal number. 

equation determines all the R^s as formal power series of the g^s upon imposing that 
Rr, = l + 0{{gk}) for aU n > 0. 

For illustration, let us consider the simple case of diagram with only bi-, tetra- and hex- 
avalent inner vertices, i.e. with gk = for all k > 4. The master equation then reads 

-Rn = 1 + 9lRn + 92Rn{Rn-l + Rn + R-n+l) 

+ gzRn{Rn-2Rn-l + R^-1 + -^n + '^Rn{Rn-l + Rn+l) (2-6) 

+ Rn-lRn+1 + Rn+1 + Rn+lRn+2) 

This relation is illustrated in Fig. 4. In the following, we will also make use of the generating 
function R for two-leg diagrams without restriction on the geodesic distance. This may be 
recovered as the limit n — > oo of i?^, hence R satisfies the relation 

^=1 + E^^('V')^' (2-7) 

k>i ^ ^ 

obtained by counting walks of length 2k — 1 with global height difference equal to — 1, 
thus having k descending steps each weighted by R. In the truncated case above, Eq.(2.7) 
reduces to 

R = l + giR + Zg2R^ + l^g^^R^ (2.8) 
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2.3. Other applications of the master equation 
The master equation may be rewritten as 

by noting that V^ ,^_i = The latter is a direct consequence of the general 

reflection relation 

Za,a-l{2k - 1) = RaZa-l,a{2k - 1) (2.10) 

obtained by first transferring the weights i?j to the ascending steps (z — 1) — > z and then 

reflecting the walks, thus interpreting them as walks from (a — 1) ^ a. Note that there 
are as many descending steps z ^ (z — 1) as ascending steps (z — 1) ^ z except for i = a 
where there is one more descent (before reflection), hence the extra factor Ra needed to 
reproduce the correct weight. 




Fig. 5: An example of labeled mobile (left) made of (inflated) polygons 
attached by their vertices into a tree-like structure. The vertex labels around 
each polygon obey the rule displayed on the right that they either decrease 
by 1 or increase weakly clockwise. 

In its alternative form (2.9), the master equation appeared in a slightly different 
enumeration problem, that of so-called labeled mobiles [12], generalizing the so-called well 
labeled trees of Ref. [13]. These mobiles are made of rigid polygons, with weight gk per 
fc-gon, glued by their vertices into a rooted tree-like object (see Fig. 5 for an example). The 
labels are subject to local rules around each /c-gon: going clockwise around a /c-gon, the 
labels must either decrease by 1 or increase weakly. This rule again is efficiently described 
by the transfer matrix Q of Eq.(2.1) by reading the labels around A;-gons from the vertex 
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at which they are suspended and translating the sequence of labels into heights forming 
a closed walk of k steps around the fc-gon. This walk is then transformed into a closed 
boundary walk of length 2k by transforming each non-negative step into a descent 

followed hj p + 1 unit ascending steps. Denoting by Rn the generating function for planted 
mobiles with root labeled n and inspecting the possible environments of the root, we have 
the relation ^ 

^ 1 - E 9kMn{k) (^-^^^ 

k>l 

where Mn{k) is the generating function for A;-gons rooted at a vertex labelled n and with 
k — 1 sub- mobiles dangling at the other vertices. In Mn{k), we have to assign a weight Ri to 
each vertex of the fc-gon labelled i, except for the root vertex. This amounts equivalently 
to assign a weight i?j to each —1 step from height i in the boundary walk, except for the 
— 1 step from the root. This leads to 

Mnik) = Zn-i,n{2k - 1) = {n\Q'"'-^\n - 1 > (2.12) 

Eqs.(2.11) and (2.12) boil down to Eq.(2.9). Again this equation must be supplemented 
with some initial conditions, for instance that i?j = for i < if we demand that the 
labels remain non-negative. 

Returning to graph enumeration problems, it was shown that mobiles with non- 
negative labels are in one to one correspondence with planar graphs with faces of even 
valences only and with a distinguished origin vertex [12,13]. The vertices of the mobile 
are in correspondence with vertices of the associated graph, and the labels represent the 
geodesic distance (on the associated graph) of these vertices from the origin vertex. In 
this framework, i?„ is now interpreted as the generating function for the graphs with a 
distinguished edge at geodesic distance less than n from the origin vertex. 

Finally, a third domain of application of the master equation comes from the study 
of spatially extended branching processes, i.e. probabilistic models for the evolution and 
spreading of a population. In this language, the index n stands for the (discrete one- 
dimensional) position of individuals and the blossom trees or mobiles are interpreted as 
the genealogical structure of families. As shown in Ref.[14], the generating function is 
then related to the probability for a population with germ at position n never to spread 
up to position 0. 

2.4- Integrahility of the master equation 

A remarkable feature of the master equation is that it may be solved exactly. More 
precisely, we will consider truncated versions of Eq.(2.5) by considering only diagrams with 

valences up to, say 2m. This amounts to set (7fc = for /c > m in which case V {Q) is an 
odd polynomial of degree 2m — 1. As explained in Ref. [8], the general solution of Eq.(2.5), 
now with arbitrary initial conditions but with the requirement that it converges at n — > oo 
takes the surprisingly simple form 

R^ = R'hl±l}hl±± (2.13) 
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where Un has a multi-soliton structure involving m — 1 integration constants Ai, . . . , Xm-i- 
Those constants may be fixed by the initial conditions. This structure is characteristic 
of (discrete) integrable systems and implies the existence of conserved quantities. More 
precisely, by inverting the system Ri = Ri{{X}) for any set z = n, n + 1, . . . , n + m — 2, 
one may in principle construct m — 1 conserved quantities for Eq.(2.5) in the form Xj = 
Aj{Rn.Rn+i, • • • , Rn+m-2) = const, independently of n. In practice, this construction is 
however difficult to implement. The remainder of the paper is devoted to the explicit 
construction of conserved quantities ab initio, i.e. without reference to the above solution 
(2.13). This construction will be carried out by use of simple combinatorial tools only. 

3. Random walks and heaps of dimers 

3.1. From random walks to heaps of dimers 

In this Section, we shall investigate a number of properties of the partition function 
of random walks Za,h{k) above. In particular, we will emphasize the connection between 
random walks and so-called heaps of dimers, leading us eventually to an inversion relation 
involving partition functions for hard dimers on a segment. All the present relations will 
be instrumental for proving the conservation statements of Section 4 below. 




-e 3^ -eg 3* 

k k 



Fig. 6: A schematic representation of the fundamental identity (3.1) for 
walks. The walks are decomposed according to either their first step (first 
line) or their last one (second line). 

A first fundamental identity is obtained by expressing the generating function for 
paths of length /c + 1 in two different ways: (i) as the concatenation of a path of length 1 
(first step) and a path of length k, or (ii) as the concatenation of a path of length k and 
a path of length 1 (last step). According to whether the first (resp. the last) step is up or 
down, we have respectively 

Za,h{k + 1) = Za+l,h{k) + RaZa-l,b{k) = Za,b-l{k) + Rb+lZa,b+lik) (3.1) 
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(a) 



stripe 



... -3-2-1 
(b) 



(C) 

Fig. 7: An example of a heap of dimers (a) and its canonical representation 
(b). The thickened dimers correspond to the right projection of the heap, 
forming in (b) a hard dimer configuration on the vertical line x = 0. By 
rotating the picture (b) and extending each dimer into a 2 x 1 rectangle, we 
obtain a piling up of these rectangles (c) whose bottom row is formed by the 
former right projection. 

as illustrated in Fig. 6. 

A second important property is that we may interpret Za,h{k) as a generating function 
for heaps of dimers, a particular case of the so-called heaps of pieces introduced in Ref. [9] . 
These heaps of dimers also occur in relation to lattice animals [10] and so-called Lorentzian 
gravity [11] and may be defined as follows. Let us decompose the plane into parallel (say 
horizontal) stripes of width 1 by drawing horizontal lines at integer vertical positions. The 
stripes are labeled by the position of their top boundary. We then place within the stripes 
a number of dimers, i.e. vertical segments of length 1 (see Fig.7-(a)). All dimers may 
freely slide within their stripe provided they do not cross a dimer within the same stripe 
or within the stripe immediately above and below. A heap of dimers is defined modulo 
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this sliding freedom and only records the relative positions of the dimers. In a canonical 
representation, we may place all the dimers at integer horizontal coordinates within their 
respective stripes with the requirement that these coordinates be negative or zero and 
maximal (see Fig.7-(b)). This is realized by pushing all dimers as much to the right as 
possible while staying at integer horizontal positions in the left half plane. This in turn 
allows upon rotating the picture by 90° clockwise to view a configuration as a piling up 
of rectangles of size 2x1 (see Fig.7-(c)) deposited on the x = line. For each heap of 
dimers, we may define its right projection as the configuration made of those dimers which 
may slide freely all the way to infinity to the right. In the canonical representation, these 
are those dimers with horizontal position 0. Clearly, the right projection of any heap of 
dimers is a configuration of hard dimers on a (vertical) line, i.e. a set of dimers occupying 
(vertical) segments [z — l,z] with the restriction that no two dimers may come in contact 
(the occupied segments must be disjoint). 

(a) 




Fig. 8: An example of a pyramid (a) and a half-pyramid (b), both with base 
[a, 6]. 

For (6 — a) odd and positive, we may consider the maximally occupied hard dimer 
configuration of the segment [a, 6]. It is made of (6 — a + l)/2 dimers occupying the 
elementary segments [a + 2z, a + 2i + 1], i = 0, . . . , (6 — a — l)/2. Any heap of dimer 
having this maximally occupied segment as right projection will be called a pyramid of 
base [a, 6] (see Fig.8-(a) for an example). We then define the generating function Ha,b{k) 
for pyramids of base [a, b] made of a total number k of dimers, with weight Ri per dimer 
in the stripe labeled z, except for the dimers of the right projection. Note that Ha,b{k) is 
non zero only if /c > (6 — a + l)/2 and that Ha^b[ip — a, + l)/2) = 1. 

We these definitions, we have the remarkable relation 



Z„,,(2A; - 1) = Ha^b{k) (3.2) 
valid for (6 — a) odd and positive, and k > 1. This relation is illustrated in Fig.9-(a). It 
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^ 

2k -1 steps 



k dimers 



Fig. 9: A schematic representation (first line) of the equahty (3.2) between 
generating functions Za^bC^k — 1) for walks from a to 6 with 2/c — 1 steps and 
Ha,b{k) for pyramids of base [a, b] with a total of k dimers. The similar repre- 
sentation (second line) for the relation (3.3) between the generating function 
Z^^{2k — 1) of positive walks and that i7^^((/c) of half-pyramids. 

follows from a general bijection between walks and heaps of pieces given in Ref. [9]. This 
bijection is recalled in detail in Appendix A below for the particular case at hand. 

In the next section, we shall also make use of the generating function Z^f^(k) for walks 
of k steps going from height a to height 6, with weight Ri per descent i —>■ [i — 1), and 
which stay at or above height a. In the following, we shall refer to these walks as positive 
walks. The generating function Z^f^{k) may be obtained from Za,b{k) by simply taking 
Ri ^ for all labels z < a. In particular, we deduce that, for (6 — a) odd and positive and 
k > 1 

ZW^k - 1) = HX,{k) (3.3) 

where H^f^{k) denotes the generating function for half-pyramids of base [a, b] with a total 
of k dimers, i.e. pyramids with dimers only in stripes with labels i > a (see Fig.8-(b) for 
an example). This relation (3.3) is illustrated in Fig.9-(b). 
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3.2. Inversion relation 

Viewing hard dimers on a line as mutually excluding and fermionic objects (with at 
most one dimer in a unit segment), heaps of dimers appear as their interacting bosonic 
counterpart in which arbitrarily many dimers may be piled up in each stripe. It is therefore 
natural to expect some boson/fermion inversion relations to hold between their respective 
generating functions. An example of such inversion relation in a "grand- canonical" ensem- 
ble, i.e. for arbitrarily many dimers, may be found in Refs.[9,ll]. In this Section, we shall 
derive another inversion formula, now in the canonical ensemble, i.e. with fixed numbers 
of dimers. 

Introducing the generating function Ila^b{k) for configurations of k hard dimers in the 
segment [a, b] and with weight Ri per dimer in the segment [^ — 1, ^], we have the remarkable 
inversion relation 

fc 

E(-l)'"'n,,«+2fc-i(/c - ^)^i+2,(2^) = (3.4) 

£=0 

for k,j > 0, with the convention that na^o-i(O) = 1. This relation may be written in a 
more compact matrix form by defining a lower-triangular, semi-infinite matrix Z(a) with 
entries 

for < j <i, and a lower-triangular, semi-infinite matrix D(a) with entries 

D(a)fc,i = {-l)''-ma,a+2k-i{k - i) (3.6) 

for < z < A; while all other entries vanish. The relation (3.4) now reads 

D(a)Z(a) = I (3.7) 

with I the (semi- infinite) identity matrix. To prove the relation (3.4), let us consider pairs 
V made of 

(i) a hard dimer configuration in the segment [a, a + 2k — 1] 

(ii) a half-pyramid with projection [a, a + 2j — 1] 

with a total number of k dimers. In the half-pyramid, each dimer in the stripe i (including 
the dimers of the right projection) receives a weight Ri. On the contrary, in the hard dimer 

configuration, each dimer in the segment [z — l,z] receives a weight —Ri. Denoting by £ 
the number of dimers in the half-pyramid, the generating function Pa{j, k) for the above 
pairs reads 

k 

Pa{j, k) = n„,,+2i-l(j) Y.{-lf-'lia,a+2k-l{k - £)H+^+,^_,{£) (3.8) 

£=0 

Note that Paij, k) is non zero only if A; > j. For each such pair V, we may absorb all dimers 
of the hard dimer configuration into a larger heap Ti{V) of k dimers by simply adding them 
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Fig. 10: An example (left) of pair P of a hard-dimer configuration in the 
segment [a, a + 2/c — 1] and a half-pyramid of base [a, a + 2j — 1] (with j < k), 
with a total number k of dimers (here k = 7). These are concatenated (right) 
into a larger heap H. The equivalence class C{H) of those pairs leading to 
the same heap Ti is constructed by considering in Ti the dimers which belong 
to the left projection of the heap but do not belong to the right projection. 
In the present example, there is exactly one such (encircled) dimer. The 
equivalence class C{H) is generated by distributing in all possible ways these 
dimers either in the hard-dimer configuration or in the half-pyramid. If we 
assign opposite weights to dimers in the hard-dimer configuration and in the 
half-pyramid, this results in a vanishing net contribution of the class C(7i) 
unless the left and right projections of Ti. are identical. This happens only 
when Ti is the maximally occupied hard-dimer configuration made of k dimers 
in the segment [a, a + 2k — 1]. 

to the left of the original half-pyramid (see Fig. 10). We can then regroup all pairs leading 
to the same larger heap TC into an equivalence class C{Ti.). Conversely, given H = TiiV), 
we may reconstruct all pairs in the class C(7i) by considering the set of dimers which belong 
to the left projection of Ti (those dimers which may be pushed all the way to infinity to 
the left) but which do not belong to the right projection. If this set contains at least one 
dimer, this dimer may be incorporated either in the hard dimer configuration of the pair 
or in its half-pyramid, contributing with opposite weights. This causes the contribution to 
Pa{j, k) of the class CiTi) to vanish. The only case where no such vanishing takes place is 
when the left projection of Ti is identical to its right projection, in which case Ti itself is a 
hard dimer configuration with k dimers on the segment [a, a -(-2/c — 1], therefore the unique 
maximally occupied hard dimer configuration, with contribution [—l)^~^Iia,a+2k-i{k) as 
the j lower dimers with label less than a + 2j — 1 belong to the half-pyramid and the [k—j) 
upper ones belong to the hard dimer configuration. We immediately deduce that 

Pa{j, k) = {-lf-ma,a+2k-l{k) Sk>j (3.9) 
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with 6k>j = 1 for A; > J and zero otherwise. Using S^j = Sk>j — Sk>j+i, we deduce 



Sk,j na,a+2fc-l(/c) = PaU, k) + PaU + 1, k) 

k 

= J](-l)^-Xa+2fc-l(/c-£) 
£=0 

X iKa+2j-ima,a+2j-lU) + <a+2, + l Wna,a+2, + l (j + 1)) 

(3.10) 

Using <a+2,-i W = <a+2,-i(2^ " 1), Ka+2j+ii^) = <a+2,+i(2^ " 1) from Eq.(3.3) 
and the relation Ila^a+2j+i (j +1) = Ra+2j+i^a,a+2j-i{j) for maximaUy occupied segments, 
we obtain 

k 

E(-l)'"'na,a+2fe-l(/c - £) ««+2,-l(2^ " 1) + ^a+2i+l<«+2i+l(2^ " 1)) = ^Kj 

£=0 

(3.11) 

which, upon using 

<a+2,-l(2^ - 1) + ^a+2, + l<«+2,+l(2^ " 1) = <a+2i(2^) (3-12) 

from Eq.(3.1), reduces to the desired inversion relation (3.4). 

Note that, as D(o) and Z(a) are lower-triangular, we may decide to truncate them 
to m X m matrices T)rn{a) and Zm(a) with indices i strictly less than m. These matrices 
clearly satisfy the inversion relation Dm,(o)Zm(o) = Im, with the mxm identity matrix. 
For illustration, taking m = 3, we have 



1 

D3(a) = I -Ra+1 1 

Ra+lRa+3 —{Ra+1 + Ra+2 + Ra+s) 1 

1 

Z3(a) = ( Ra+1 1 

Ra+1 {Ra+1 + Ra+2) {Ra+1 + Ra+2 + Ra+3) 1 



(3.13) 



which are inverse of one-another. 



4. Conserved quantities 

4.1. Definition of the conserved quantities 

We now present compact expressions for a particular "basis" of conserved quantities 
of the master equation (2.5). More precisely, we shall define below a set of quantities 
{^2i{n)} for i > 1 and n > satisfying r2i{n)— const, independently of n. Of course, any 
functions of the r2i's are conserved as well and the precise choice of basis below has been 
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done in regard of its combinatorial nature. Other, alternative choices will be discussed in 
Sects. 6.1 and 6.2 below. 

Using the generating function Z^^{k) of previous section for positive walks and the 
grand-canonical generating function for walks of odd length, we define the quantity 
T2^{n) by 



^n-l,n-l(2^)ro(n) 



^n-l,n-l+2j 

3 = 1 



+2j-l,n-2 



(4.1) 



for i > 1 and n > 0, with ro(n) given by 



ro(n) = Rn-l - V^_^ 2 + ^n,0 



(4.2) 



for n > 0. The quantity ro(n) itself is not stricto sensu a conserved quantity. However, the 
master equation precisely ensures that Voin) = 1 for n > 1 while ro(0) = 1 by definition 
(since R-i = V^i _2 =0), hence the conservation of Fq is a tautology. This in turn allows 
to simplify the expression (4.1) for r2i(n) by substituting ro(n) = 1. The proof of the 
conservation of r2i(n) however is made much simpler by keeping as such the slightly more 
involved definition (4.1). Introducing the vectors T{n) and V'{n) with components 



r(n)i = T2i{n) and V'{n)j = {Rn-\ + K,o)5j,Q - K+2j-i,n-2 



(4.3) 



for > 0, the above relations (4.1) and (4.2) read simply 



f(n) = Z(n - 1) V'{n) 



(4.4) 



with Z(a) defined as in (3.5). It will prove useful to invert this relation. This is readily 
performed by use of the inversion relation (3.4) or (3.7), with the result 



V'(n) = D(n-l)f(n) 



(4.5) 



with D(a) defined in (3.6), namely 



K+2,-l,n-2 = Y.^-iy-^^n-l,n+2j-2{i)T2j-2i{n) + o(i?n-l + 5n,o) (4.6) 
i=0 



with again the convention that n^_i „_2(0) = 1. 
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4-2. Proof of the conservation 



We are now ready to prove that the quantities T2i{n) of Eq.(4.1) are conserved quan- 
tities of the master equation (2.5). To this end, we shaU now use the fundamental equation 
for V 



K'+l,6 + = <6-l + i?6+l<6+l (4-7) 



directly inherited from the fundamental equation (3.1) for Za,h{k)- Taking a = n -\-2j — 1 
and 6 = n — 1, we get the identity 



^(n+l)+2i-l,(n+l)-2 + -Rn+2i-l V('„+i)+20-l)-l,(n+l)-2 

= K+2j-l,n-2 + -Rn^(n+2)+2(j-l)-l,(n+2)-2 



For J > 1 and n > 0, we may substitute Eq.(4.6) into Eq.(4.8), leading to 



( - T2j{n + 1) + Y,{-lr-^Iin,in+l)+2i-2{^)T2J-2^{n + 1)) 
1=1 

+Rn+23-l{^^{-'^y~^'^n,{n+l)+2{3-l)-2{i)^2{3-l)-2i{n+ 1) + 5j_i,o-Rn) 
i=0 
3 

= ( - r^jin) + Y,{-iy-'Un-l,n+2J-2ii)r2j-2i{n)) 
i=l 

+-Rn( y^(-l)'~'^n^+l,(n+2)+20-l)-2(^)r2(j-l)-2i ("■ + 2) + 5j_i,o-Rn+l) 



i=0 



where we have used Iln-i^n+2j-2{0) — n^^„_|_2j-i(0) = 1 for all j > (unique empty dimer 
configuration). Noting that the boundary terms on both sides cancel, we now change 
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i ^ (z — 1) in the last sum of each side and rearrange the factors of F, leading us to 



T2jin + 1) -T2jin) 
j 



n,n+2j — 1 

2=1 

j 

- y^(-l)'~'^{n„_i,n+2j-2(0 - ^nnn+l,n+2i-2(« - l)}r2j-2i(n) 
i=l 

i 

- J](-l)*i^,^„+l,„+2,-2(^ - 1) {T2J-2^{n + 2) - r2,_2,(n)) 
1=1 

j 

Y,{-^y~^^n,n+2j-2{l) {V2J-2^{n + 1) - r2,-22(n)) 

:1 

3 

Y,{-iyR^Iir,+l,n+2j-2{i " 1) {V2j-2r{n + 2) - r2,_2z(n)) 



2=1 



i=l 



(4.10) 
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Fig. 11: A schematic representation of the fundamental identity (4.12). A 
configuration of hard dimers on the segment [a, h] may be viewed either as 
a configuration on the segment [a — 1 , 6] where the unit segment [a — 1 , a] is 
empty (hence the subtraction) or as a configuration on the segment [a, 6 + 1] 
where the unit segment [6, 6+1] is empty. 

In the last equality, we have used the property 

nn,n+2j-l(0— -Rn+2j-inn,n+2j-3(^ ~ 1) 

= n„-l,n+2j-2(^) — RrJ^n+l,n+2j-2{'i' ~ 1) = n„,n+2j-2(^) 



(4.11) 
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which is a particular instance (a = n, b = n + 2j — 2) of the fundamental identity 



Ua,b{i) = Ila-i,b{i) - Ralia+l,b{i - 1) = - Rb+llla,b-l{i ' 1) (4.12) 

satisfied by generating functions for hard dimers. This identity is illustrated in Fig. 11 and 
may be viewed as the hard dimer counterpart of Eq.(3.1) for random walks. 
Prom Eq.(4.10), we immediately deduce by induction on j that 

r2,(n+l) = r2,(n) (4.13) 

for all j > I and n > provided that ro(?i + 1) = ro(n) for all n > 0. As already 
mentioned, this last requirement is precisely guaranteed by the master equation (2.5). The 
T2iS therefore form a set of conserved quantities for this equation. 

For illustration, we list below the first two conserved quantities for the truncated case 
of graphs with valences up to 2m = 6, corresponding to the truncated master equation 
(2.6): 



with 



r2(n) = i?^-v;:+i,„_2 

r4(n) = Rn{Rn + Rn+l) - (Rfi + Rn+1 + -Rn+2)K+l,n-2 " K+3,n-2 

V^+l,n-2 — Rn+lRnRn-1 ((72 + g^{Rn+2 + -Rn+l + -Rn + Rn-1 + Rn-2)) 
yn+Z,n-2 — 93Rn+3Rn+2Rn+lRnRn-l 

where we explicitly substituted To{n) = 1 into Eq.(4.1). 



(4.14) 



(4.15) 



5. Graph interpretation 

5.1. Conserved quantities as multi-point correlation functions 

The constant value of the conserved quantities has a nice combinatorial interpreta- 
tion as a multi-point correlation function. More precisely, let us define by G2i{{gk}) the 
so-called (disconnected) 2i-point function, i.e. the generating function for possibly dis- 
connected 2i-leg diagrams, namely graphs with inner vertices of even valences (weighted 
gk per 2/s-valent vertex) and with 2i legs, i.e. univalent vertices, adjacent to the same 
(external) face. These legs are distinguished and labeled 1, 2, . . . , 2z counterclockwise. We 
have the identification 

^2^(n) = G2i (5.1) 

for all ?i > and all z > with the convention that Gq = 1. To prove this, we simply 
evaluate r2i(n) at n = 0. Noting that V2j_i _2 = as it is proportional to = 0, we 
are left with 

r2i(0) = ^^l,_l(2^) = Z_i,_i(2z) = Zo,_i(2z - 1) (5.2) 
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(a) 



(b) 



(c) 



Fig. 12: Bijection between (a) two-leg diagrams witti botii legs in the same 
face and whose outcoming leg is adjacent to a 2i-valent vertex and (c) 2i-leg 
diagrams. In the intermediate step (b), we have glued the two legs into a 
rooted edge and cut out the 2z-valent vertex. The legs are labeled counter- 
clockwise by 1, 2, . . . 2i with the former incoming leg labeled 1. Note that the 
2z-leg diagrams need not be connected in general. 

where we first note that R-i = automatically selects positive walks and we then remove 
the first (ascending) step with weight 1. 

As explained in Sect. 2. 2, the quantity Zo,-i(2i — 1) is the generating function for 
blossom trees with a 2z-valent root vertex and a contour walk with depth. Returning 
to the graph interpretation, it is identified via the bijection of Sect. 2.1 as the generating 
function for two-leg diagrams with an outcoming leg attached to a 2i-valent vertex and 
adjacent to the same (external) face as the incoming leg. The 2i-valent vertex receives a 
weight 1 as opposed to all other inner vertices weighted by the usual factors gk if they are 
2A;-valent. Such two- leg diagrams are in bijection with the 2z-leg diagrams defined above. 
Indeed, starting from these two-leg diagrams and gluing their two legs into a marked edge, 
we simply erase the unweighted 2i-valent vertex and obtain the desired 2z-leg diagrams 
with a marked leg which receives the label 1 (see Fig. 12 for an example). Note that 2i-leg 
diagrams need not be connected except for i = 1 and may split into connected 2j-leg 
diagrams with some j's summing to i. Eq.(5.1) follows. 

An outcome of this result is that we obtain a host of explicit expressions for G2i in 
terms of the generating functions Ri of blossom trees corresponding to various choices of 
n in Eq.(5.1). The first expression is nothing but G2i — Zo,-i(2i — 1) corresponding to 
n = 0. For illustration, for z = 1,2 and 3, we have 



Another particularly simple choice consists in using Eq.(5.1) in the limit n — > cxd. This 
gives explicit formulas for G2i in terms of the function R solution of Eq.(2.7) (recall that 



G2 — Rq 



G4 — Ro{Ro + Ri) 

Gq = Ro(Rq + 2RqRi -\- Ri + R1R2) 



(5.3) 
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R is simply the generating function for two-leg diagrams with arbitrary geodesic distance 
between the legs). Letting — > i? for all a, we first find that 



(5.4) 



k>j+ 

which finally yields from Eq.(4.1) 



mi 



fc>2 j=l 



i- jj V - j - J \ k + j 



These expressions are determined equivalently as the solutions of the inverse relation (4.6) 
at large n, which reads 

where (^"'~*) is the number of hard dimer configurations of i dimers on a segment of length 
2j-l 

For illustration in the truncated case with valences up to 2m = 6, we find from Eq.(5.5) 
that 

G2 = R- g2R^ - 5g3R'^ 

G4 = 2i?2 - 3g2R'^ - IQgsR^ (5.7) 

Gq = 5R^ - 9g2R^ - ^OgsR^ 

where R is now determined by Eq.(2.8). 

The explicit expressions (5.5) or (5.6) are to be compared with other known expressions 
for multi-point correlation functions. The 2z-point function above may indeed be computed 
alternatively cither by use of the planar limit of the one matrix model or by the so-called 
loop equations. The planar solution of the one-matrix integral may be obtained via saddle 
point techniques and reads in the so-called one-cut case: 

i2^2i-2e(^^]R' = lim {n\Q^'{Q - V'{Q))\n - 1) 

' \ t J n— >oo 

■ r-. u 5.8 

'^ty-j:.^- E a;; 

^ k>2 j^l ^ J / \ 'J 
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This expression is readily equivalent to Eq.(5.5) by simply noting that 



E 



2A / / 2i - 2£ \ / 2i-2i \\ _ /'2i + 1 
e) [[i-i-jj ~ - [i-j 



(5.9) 



obtained by cutting any walk of length 2i + 1 from, say to height 2j + 1 at the level of 
its last — > 1 step, resulting into a first walk of length, say 2£ from to and a positive 
walk of length 2{i — £) from height 1 to height 2j + 1. 

On the other hand, the loop equations simply express the (disconnected) 2i + 2-point 
function as a sum of either a 2z + 2/c-point function if the first leg is connected to a 2A;-valent 
vertex or a product of two lower order correlations if the first leg is connected directly to 
the {2j + 2)-th leg, namely: 

i 

G2i+2 — gkG2i+2k + G2i-2jG2j (5.10) 
k>l j=0 

That Eq.(5.5) solves this loop equation may be seen as a necessary consistency of the 
planar limit of the one-matrix model in the one-cut case which implies that Eqs.(5.10) and 
(5.8) are compatible. 

5.2. Relations between conserved quantities for bounded valences 

In the truncated case just above of valences up to 2m = 6, Eqs.(5.7) imply the relation 
giG2+g2G4+g3GQ + l = G2 provided Eq.(2.8) is satisfied. This may be seen as a particular 
case of the loop equation (5.10) for z = above in its truncated form. More generally, the 
truncated loop equations for graphs with valences up to 2m allow for expressing all G2i 
for i > m as polynomials of the first m — 1 values G2, G4, . . . , G2(m-i)- This gives a set of 
polynomial relations between the values of the conserved quantities T2i{n). 

Using the results of previous Section, this interdependence may be rephrased into 
linear relations between the G2i's by writing Eq.(5.6) for j > m in which case the l.h.s. 
vanishes. This results in 

0- J](-l)^-^f^^p)i?^G2,-2i (5.11) 



i=0 



where R is determined by the polynomial truncation of Eq.(2.7) with gk = Q for k > m. 
For illustration, for m = 3, we have the first two relations 

Ge = 5RG4 - m^G2 + R^ 

(5.12) 

Gs = TRGq - I^R^Ga + im^G2 - R^ 
and similar relations for G2i's with higher indices. 
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Fig. 13: Any blossom tree (a) whose root is encircled by z > 1 edges in the 
closing process of Sect. 2.1 may be rerooted at the first excess bud clockwise 
from the root, while this original root is replaced by a leaf (b). This results 
in a rooted tree with a, say 2/c-valent root vertex adjacent to /c — 2 buds and 
k + 1 blossom subtrees (with necessarily k > 2). The root of this new tree 
is encircled in the closing process by i — 1 edges. In (c) and (d) we focus on 
the rearrangement of bud-leaf pairs in a general case. The excess buds in (d) 
are those of (c) except the first one, now promoted to root. This leaves three 
unmatched leaves. 

5.3. Combinatorial interpretation o/r2(n) 

The conservation of r2(n) has a nice combinatorial explanation in the language of 
blossom trees. Indeed, assuming Toin) = 1 for all n > 0, the equality T2{n) = T2{0) for 
all n > 1, with r2(n) defined by Eq.(4.1), may be rewritten as 

(i?„-i?o) = <+i,n-2 (5-13) 

The quantity Rn — Ro is the generating function for blossom trees with contour walk of 
depth i with < z < n. In other words, in the closing process of these trees into two-leg 
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diagrams, their root is encircled by at least one and at most n edges separating it from 
the external face. Picking the bud from which the deepest encircling edge originates, we 
may reroot the tree at this bud and replace the original root by a leaf (see Fig. 13). The 
resulting object is a tree satisfying (Bl), (B2) and (B3) except for the root vertex which 
now has exactly {k — 2) buds if it is 2A;-valent (note that by construction, this vertex has 
valence at least 4 as it originally carried a bud). In particular, all the proper subtrees 
of this tree are ordinary blossom trees. Finally, the contour walk of this new tree (now 
stepping from height to height +3) has clearly depth (i — 1) < (n — 1). The desired 
generating function may again be obtained by reading the sequence of buds and subtrees 
counterclockwise around the root vertex and reads Zn+i^n-2{2k — 1) if the new root vertex 
is 2fc-valent. Summing over all possible valences, we get the generating function l^+i^^_2- 
As the above re-rooting procedure clearly establishes a bijection between the two types of 
trees at hand, this provides a bijective proof of the equality (5.13). As this equality is valid 
for all n > 1, this in turn proves the conservation of r2(n). 

It would be nice to have similar bijective proofs for the conservation of T2i{n) for z > 2 
as well. For 2i = 4, it is easily seen that, upon using To{n) = 1 and r2(n-|-2) = r2(n) = Rq, 
the equality T4^{n) = ^^{O) may be rewritten as 

Ro{Rn+l - Rl) = Vn+3,n-2 + ^n-|-3,n^n-|-l,n-2 (5-14) 

The l.h.s. of this identity is nothing but the generating function for pairs made of a blossom 
tree with contour walk of depth and a blossom tree with contour walk of depth i with 
1 < i < n + 1. The r.h.s. may be interpreted as the generating function for the collection 
of two types of objects 

(1) rooted trees satisfying (Bl), (B2) and (B3) except at the root vertex which carries 
exactly {k — 3) buds if it is 2A;-valent (with necessarily k > 3). The contour walk of 
this tree moreover has depth at most n — 1; 

(2) pairs of rooted trees satisfying (Bl), (B2) and (B3) except at their root vertex which 
carries exactly {k — 2) buds if it is 2/c-valent (with necessarily k > 2). The first of 
these trees has a contour walk of depth at most n + 1 and the other a contour walk 
of depth at most n — 1. 

This suggests the existence of a bijection between the two collections of objects above. If 
such a bijection could be exhibited, it would provide an alternative bijective proof of the 
conservation of r4(n). 

For larger i, the conservation of r2i{n) yields more involved relations of the type 
above which may suggest higher order bijections as well. For instance for 2i — 6, we get 
the relation 

RoRl{Rn+2 — R2) — Ro{Rn+l — Rl){Rn+3 — Rl) = 

Ki+5,n-2 + Ki+5,nKi+l,n-2 + ^n+5,n+2^n+3,n-2 + Ki+5,n+2Ki+3,nKi+l,n-2 

(5.15) 
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while for 2z = 8, we have 

RoRlR2{Rn+3 — R3) — RqRi (j^Rn+l — -Rl)(-Rn+4 — R2) + {Rn+2 — -R2)(-Rn+4 — Rl) 
+ (-Rn+2 — R2)iRn+5 — -^2)^ + Ro(j,Rn+l — -Rl)(-Rn+3 " -Rl)(-Rn+5 " Rl)j 

— Ki+7,n-2 + Ki+7,nKi+l,n-2 + ,n+'^n+b,n-2 

+ ^n+7,n+2^n+3,n^n+l,n-2 + ^n+7,n+4^n+5,n^n+l,n-2 + ^n+7,n+4Ki+5,n+2K(+3,n-2 

_i_ T/' T/' -I- T/' V' V' V' 

^ ^n+7,n+2^n+3,n-2 ' ^ n+7,n+4:^ n+5,n+2^ n+3,n^ n+l,n-2 

(5.16) 

These relations still await a bijective proof. 
6. Discussion 

6.1. Alternative expressions for the conserved quantities 

As already mentioned, the definition (4.1) corresponds to a particular choice of con- 
served quantities involving positive walks. This choice however does not respect the "time 
reversal" symmetry of the master equation, namely under Rn-j Rn+j for ah j (see 
Eq.(2.6) for illustration). In this Section, we present another choice of conserved quanti- 
ties that respect this symmetry. More precisely, we define 



^2i{n) = {Zn,n-l{'2i - 1) - Rn-lRn+lZn-2,n+l{'2i " i)}'^o{n) 

i 

- {-^n-j,n+i-l(2i - 1) - Rn-j-lRn+j+lZn-j-2,n+j + l{'^'i - 1)} V^n+i,n-i-l 



(6.1) 

for z > 1 and n > 0, with ro(n) given by 

fo(n)=i?.-y;_i (6.2) 

for n > 0. Note that in both Z's and 's, the indices lead to symmetric expressions. 
By techniques similar to those of Sect. 4.1, one can invert Eqs. (6.1) and (6.2) into 

j 

VUj,n-i-l = E(-l)'"'nn-,-,n+i-l(i)f 2i-2i(n) + 5,, oi?n (6.3) 

with the convention n„^„_i(0) = 1. The passage from Eq.(6.1) to Eq.(6.3) makes use of 
an inversion relation similar to Eq.(3.4). This relation is explicited in Appendix B. 

The conservation of T2i{n) may be derived along the same lines as in Sect. 4. 2. Using 
again Eq.(4.7) now with a = n + j and b = n — j, we have the relation 

^(n+l)+i,(n+l)-j-l + Rn+jV^+(j-i),n-ij-l)-l 

= ^n+j,n-j-l + Rn-j+lV[n+l) + {j-l),in+l)-ij-l)-l 
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For j >1 and n > 0, substituting Eq.(6.3) into this relation yields 

j 

i=l 

J-1 

+Rn+j{^{-iy~^'n.n-{j-l),n+U-l)-li^'^2{j-l)-2i{n) + 5(j_i),oi?n) 
1=0 

J 

= ( - f 2,(n) + ^(-l)^-in„_,- „+,_i(z)f 2,_2i(n)) 



J-1 



^n(^+i)_(j_i)^(„+i) + (j_i)_i(*)f2(j_i)_2i(n+ 1) +%_l),oi?n+l) 



i=0 



(6.5) 

As the boundary terms on both sides cancel, upon setting i — > {i — 1) in the last sum of 
each side, we may rewrite this identity as 

f2,(n+l)-f2,(n) 

j 

= ^{-'^y~^ {'nn-j+l,n+j{^) + Rn-j+l'^n-j+2,n+j-l{i " 1)} f 2j-2i + 1) 
1=1 

J 



- ^ {nn-j,n+j-l(*) + -Rn+jn„_j + i,„+j_2(i - l)}T2j-2i{n) 

i=l 

= {Sn-,+l,n+i-l«} (f 2,-2i(n + 1) - f 2,-2i(n)) 



i=l 

(6.6) 

where we have defined 

In the last line of Eq.(6.6) we have used the relation Ila^b+iii) + Ra^a+i,b{'>' ~ 1) = 
^a-i,bii) + Rb+i^a,b-i{i ~ 1) = ^a,b{'i-) which is a direct consequence of Eq.(4.12). We 
immediately deduce from Eq.(6.6) that all T2ji'n) for j > 1 are conserved provided ro('^) 
is independent of n. This is the case when the master equation (2.5) is satisfied, as it reads 
precisely ro('^) — 1 ^or all n > 0. This completes the proof of conservation of r2i- 

For illustration, we list below the two independent conserved quantities for the trun- 
cated case with up to 2m = 6-valent graph corresponding to the truncated master equation 
(2.6): 

f2(n) = Rn- K+l,n-2 

r4(n) = Rn{Rn-l + Rn + Rn+l) — Rn-lRn+1 (6-8) 
- (Rn-1 + Rn + Rn+l)V^_^.i n-2 " K+2,n-3 
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with 

V^,ln_2 = Rn+lRnRn-l{g2 + gd,{Rn+2 + Rn+1 + -Rn + Rn-1 + Rn-2)) 

(6.9) 

V„+2,n-3 — 93Rn+2Rn+lRnRn-lRn-2 

where we exphcitly substituted ro(^) = 1 into Eq.(6.1). 

When compared with Eq.(4.14), we see that T2{n) = ^2{n) and r4(n) = r4(n — 1) + 
{Rn+i + Rn + Rn-i)(^2{'n') — T2{n — l)). This shows that the new set of conserved quantities 
is not independent from that of Section 4, as expected. 

This interdependence can be made even more explicit by noting that, for n = 0, we 
again have r2i(0) = Zo,_i(2i — 1) = r2i(0). When Eq.(2.5) is satisfied, the constant value 
of r2i('^) is therefore equal to G2i- This is also apparent from the large n limit of Eq.(6.3) 
which is identical to the large n limit of Eq.(4.6). As a last remark, wc note that the 
conservation of the F's is granted by that of the F's as we may then replace r2j-2i{n) in 
the r.h.s. of Eq.(6.3) by a shifted value r2j-2i{n — j ' + 1) which, upon a global shift of 
n —>■ n + j — 1, show that the r2i's also obey Eq.(4.6) for j > I. Together with the master 
equation which guarantees that To{n) = ro(n) = 1, this implies that, when the F's are 
conserved, we necessarily have T2i{n) = T2i{n) for all n, hence these are conserved as well 
and take the same constant values. 



6.2. Compacted conserved quantities 

Let us again restrict ourselves to the truncated case of valences up to 2m and con- 
centrate on the (m — 1) first conserved quantities r2i(^^) for i = 1, . . . , m — 1. From the 
definition (6.1), we note that the largest index carried by an Ra comes from the contribu- 
tions 5frra-^n+j,n-j-i(2m — 1) in the term V^_^_j^^_j_i appearing in the r.h.s.. This index 
is equal toa = n-|-m — 1 and reached for each j by a unique path of length 2m — 1 

starting with m — j — 1 up steps followed by m + j down steps. This term reads precisely 
j 

~9m n Rn-e- Similarly, the largest index carried by an Ra in Fo(n) comes from the 

£=l-m 



contribution — 5f^Z„^„_i(2m — 1) in the term —V^ and reads —gm Yl Rn-£- 

£=l-m 

We have the remarkable identity 

-^n-l,n-l(2i) = {Zn,n-l{^'i " 1) " -Rn-l-Rn+1 -^n-2,n+l (2i — 1)} 

+ {-^n-j,n+j-l(2* — 1) — Rn-j-lRn+j+lZn-j-2,n-\-j+l {2i-l)}l[Rn-£ 

j=l £=1 

(6.10) 

which may be inverted into 

J j 

- H Rn-£ = Y.^-iy~^Iln-^,n+j-l{i)Zn-l,n-l{2j - 2i) (6.11) 
£=1 i=0 
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These identities may be proved with arguments similar to that of Sect. 3. 2 and foUow from 
general inversion formulas listed in Appendix B. We deduce that in the combination 



Q2i{n) = t^iin) + (l - fo(n)) Z„_i,„_i(2z) 



(6.12) 



all terms containing Rn+m-i are cancelled. Note that the definition (6.12) does not involve 
the particular value of m and therefore gives a universal recipe for compacting the conserved 
quantities, irrespectively of the precise (finite) order of truncation. The quantities &2i{'n) 
form clearly an alternative set of conserved quantities that take the same constant values 
G2i, but that involve one less Ra than the original or r2i. 
For illustration, for m = 3, we have 



Q2{'n) = Rn + Rn-l — RnRn-l ^ — 91 — g2{Rn + Rn-l) 

~ 93 {Rn+lRn + {Rn + Rn-l)^ + Rn-lRn-2 — Rn+lRn-2) 
04(n) = Rn+lRn + Rn-lRn-2 + (Rn + Rn-l)^ 



Rn Rn—1 



(1 — gi){Rn+l + Rn + Rn-l + Rn-2) 



— 92{Rn+l + Rn + Rn-l){Rn + Rn-l + Rn-2) 

— gsiRn+l + Rn + Rn-l + Rn-2) {Rn+lRn + (Rn + Rn-l)^ + Rn-lRn-2) 



The value of O2 (n) at (/i = (73 = was already obtained in Ref. [14] . 
6.3. Application to tetra- and hexavalent graphs 



(6.13) 



In the truncated case of graphs with valences up to 2m, the truncated master equation 
is a recursion relation on n allowing for the recursive determination of all R^s with n > m— 
1 from the knowledge of Rq, . . . , Rm-2- These initial values are completely determined 
by using the m — 1 first conserved quantities in the form r2i(0) = r2i(oo) = G2i (or similar 
equalities for or ©2i) with z=l,...,m— 1. For instance, in the case m = 3 of graphs 
with bi-, tetra- and hexavalent vertices only, this yields, by use of the first two lines of 
(5.3) and (5.7) 



Ro = R~ 92R^ ~ ^gsR"" 

R{1 - QgsR^ - glR'' - 25giR^ - g2{R^ + lOgsR^)) 



Ri = 



(6.14) 



1 - g2R^ - hg^R^ 



with R solution of Eq.(2.8). This result is in agreement with the expressions obtained in 
Ref.[8] via the exact solution (2.13). 

Another possible application concerns local environments within graphs. For instance, 
we have access to the distribution of the number of faces adjacent to the external face in 
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rooted planar graphs, namely graphs with a marked oriented edge, with the external face 
adjacent to this rooted edge and lying on its right. These rooted graphs are in bijection 
with two-leg diagrams with both legs in the same (external) face, as readily seen by gluing 
the two legs counterclockwise around the diagram into a rooted edge oriented from the 
outcoming leg to the incoming one. The rooted graphs are therefore enumerated by Rq. 
We may then assign an extra weight x to each face adjacent to the external one, i.e. a 
weight whenever the external face has p adjacent faces. The corresponding generating 
function, still denoted Rq = Ro{x), is determined by the master equation (2.5) for all n > 1 
and by a modified equation at n = 0, namely 

i?o = a; + Fo'-i (6-15) 

In particular, all conserved quantities are valid for n > 1 and their constant value remains 
unchanged, given by G2i- 

In the truncated case of valences up to 2m, we obtain a closed algebraic system 
involving Rq, Ri, . . . , Rm-i by supplementing Eq.(6.15) by the m — 1 first conservation 
laws ©2i(l) = for i = 1, . . .m — 1. Upon eliminating all -R„'s with n > 0, we are left 
with a single algebraic equation determining Ro{x). Let us now illustrate this in the cases 
of pure tetravalent and pure hexavalent graphs. 




Fig. 14: The twelve rooted tetravalent planar graphs with up to 2 vertices, 
hence with weights p — 0,1,2. Assigning a weight x to each face adjacent 
to the external one, we reproduce the first three terms of the expansion (6.17) 
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For pure tetravalent graphs, we obtain the equation 



x{x - 1) + {1 - X + g2R{R - 2) - 2glR^)Ro + xg2Rl = (6.16) 

with R solution of i? = 1 + 3g2R^. This leads to the expansion 

Rq = x + g2{x + x^) + c/|(3a; + Ax'^ + 2x^) + gi{Ux + 20x^ + 15x^ + 5a;^) + . . . (6.17) 

The corresponding rooted graphs up to order gf| are displayed in Fig. 14 for illustration. 
From the knowledge of Ro{x), we may extract the distribution of the number of faces 
adjacent to the external face in large tetravalent rooted graphs by considering the limit 

Aix) ^ Yl ^"nP) = lim ^TT^f^ (6.18) 

^ — ' iV^OO itnl ill JV 

Here V{p) denotes the probability for the external face of infinitely large tetravalent rooted 
graphs to have p adjacent faces. The large N limits above are governed by the approach 
to the critical point g2 — g2 = 1/12 where R and Rq become singular. Expanding in 
powers of e = -\/(^^"^~"^2)7^' above ratio giving A(a;) is obtained as the ratio of the 
coefficients of in the expansions of Ro{x) and -Ro(l) respectively. Using Eq.(6.16), we 
get the following algebraic equation for A (a;) 

27x{x - 1) + 4(4 - 3a;)^A(l - xA) = (6.19) 

hence 

A(x) = , _ 

2x\ (4-3a;)i 

and 

For pure hexavalent graphs, we obtain the equation 

x{x-lf-{x-l){x-l + 2gsR^i3-2R) + 24glR^)Ro+Rxg3{R-2-5g3R^)Rl+x^g3Rl, = 

(6.22) 

with R solution of i? = 1 + 10^3-R^. This now leads to the expansion 

Ro^x + g3{2x + 2x^ + x^) + gl{28x + 32x^ + 25x^ + 12x^ + Sx^) + . . . (6.23) 

From Ro{x), wc may again extract the distribution of the number of faces adjacent to the 
external face in large hexavalent rooted graphs by considering the same limiting ratio as 
in (6.18) with g2 (73 and where V{p) now denotes the probability for the external face 
of infinitely large hexavalent rooted graphs to have p adjacent faces. The corresponding 
large N limits are now governed by the approach to the critical point gs ^3 = 2/135 



A(^) = 7^ fl - ) (6-20) 



32 



where R and Rq become singular. Performing a similar expansion as before, we arrive at 
the following algebraic equation for A (a;) 

125x{x - if - 9{x - l)(125a;"'^ + A75x^ + 200x - 96) A - 27x{x + A) (7 - 5xfA^{l - xA) = 

(6.24) 

from which we deduce the probabilities 

^/-.N 125 1625 865625 

V(l) = , V(2) = , P(3) = (6.25) 

^ ^ 864' ^ ^ 10368' ^ ^ 5971968 ^ ^ 

for the external face to have 1, 2 or 3 neighboring faces. 



7. Conclusion 

In this paper, we have explicitly constructed various sets of conserved quantities for 
the (integrable) master equation which determines the generating function R^ for two-leg 
diagrams with vertices of even valences. These conserved quantities are constructed in 
terms of generating functions for random walks, whose relations with heaps of dimers and 
hard dimers were instrumental to grant the conservation. 

The precise properties that we used, namely some boson/fermion correspondences and 
their associated inversion relations, as well as some fundamental relations satisfied by ei- 
ther walks or hard dimers, are very general and presumably can be adapted to other cases 
of interest. Indeed, in the planar graph enumeration framework, it was already recognized 
that many other classes of (decorated) graphs have generating functions which satisfy inte- 
grable systems of algebraic master equations. These include graphs with arbitrary even or 
odd valences, as well as the so-called p-constcUations considered in Ref. [15]. In both cases, 
the two-leg diagram generating functions display a form similar to that of Eq.(2.13) [8,16]. 
In a particular case of 3-constellations, namely that of graphs with bicolored trivalent ver- 
tices (so-called bicubic maps), a conserved quantity was already identified in Ref. [14]. We 
suspect that the above methods can be extended to the even more general case of planar 
bipartite graphs, also enumerated (without geodesic distance) by two-matrix models. 

Beside the dependence on the geodesic distance, we know from matrix-model analysis 
that another integrable structure underlies the topological expansion for the generating 
functions of graphs with arbitrary genus. It would be desirable to relate the two directions 
of integrability. A general theory of discrete integrable equations was built in Ref. [17] with 
(Hirota bilinear) master equations involving more than one integer index n (see also [18] 
for the theory of discrete Painlcvc equations). This suggests to look for a larger integrable 
structure that would include both the dependence on (possibly several) geodesic distances 
as well as that on topology. 

As explained in Sect. 5. 3, the conservation properties have a strong "bijective flavor" in 
the language of blossom trees. A bijective explanation for the conservation would be very 
instructive and seems to be within reach in view if the relations (5.15) and (5.16). If such 
a proof exists, it must involve only generic properties of blossom trees, or equivalently of 
planar graphs. This could allow to trace back the somewhat mysterious generic appearance 



33 



of conserved quantities in the context of planar graph enumeration problems to these 
generic properties. 

Finally, we may hope that the conservation properties have a natural interpretation 
of their own both in the language of labeled mobiles and in that of spatially extended 
branching processes. This is yet to be found. Note in this context the recent appearance 
of yet another integrable master equation governing the so-called "embedded binary trees" 
studied in Ref.[19]. 

To conclude, the general case of graphs with vertices of even valence is known to give 
access to multicritical transitions points of 2D Quantum Gravity by proper fine tuning of 
the parameters gi. As we have now a complete picture of the conserved quantities, we may 
use them to describe these multicritical points. In the scaling limit of large graphs, the 
master equation is known to reduce at such points to an integrable differential equation 
related to the KdV hierarchy [8]. Our conserved quantities provide explicit discrete coun- 
terparts of the integrals of motion for these equations. These integrals are indeed recovered 
from our expressions in the scaling limit. 
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Appendix A. Relation between heaps and walks 

We consider the generating functions Za,b{'2k — 1) for random walks of 2/c — 1 steps 
from height a to height b with weight Ri per descent i {i — 1), and its truncated version 
Z^f^{2k — 1) restricted to "positive" walks with heights larger or equal to a. Similarly, we 
consider the generating function Ha,b{k) for pyramids of k dimers with base [a, b] (with 
(b — a) odd and positive) with weight Ri per dimer in the stripe i except for those in the 
right projection. We also consider its truncated half-pyramid version H^j^{k) with dimers 
only in stripes i > a. In both cases, the truncation simply amounts to taking Ri ^ for 
all i < a. 

We wish to prove that 

Za,ti2k - 1) = Ha,bik), Z+,i2k - 1) = H+,ik) (A.l) 

for all /c > 1 and (6 — a) odd and positive. 

These relations are a consequence of the following bijection, borrowed from Ref.[9], 
between on the one hand random walks of 2A; — 1 steps from height a to height b > a and 
on the other hand pyramids with base [a, b] and a total of k dimers. Starting from a walk, 
as illustrated in Fig.l5-(a), viewed as the time (horizontal) evolution of a walker on the 
integer (vertical) line, the walker acts as the Little Thumb by dropping small pebbles at 
each passage in a unit segment unless a pebble is already present in which case he picks 
it up. As shown in Fig. 15- (b), putting dimers at the times/locations of the pick-up events 
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Fig. 15: Bijection between walks from a to b > a with 2k — 1 steps and 
pyramids of base [a, 6] witli k dimers, as explained in the text. In (a), the 
times/locations of pebble droppings are indicated by filled black ellipses. Each 
pebble stays in place during some time (indicated by a black horizontal seg- 
ment) until it is possibly picked up. In (b), we put a dimer at each pick-up 
time/location. This configuration is completed into a pyramid of base [a,b] 
in (c). The inverse mapping makes use of the splitting of dimers into up- or 
down-pointing triangles (d), whose sequence along each line determines the 
steps of the walk. 

results in a heap of dimers. Moreover, if the walker goes from a to b > a, there are b — a 
pebbles left behind. By completing the heap on the right by the maximally occupied hard 
dimer configuration on the segment [a, 6], we obtain a pyramid of base [a, 6] and with a 
total of k dimers if the walk has 2k — 1 steps, as illustrated in Fig. 15- (c). This establishes 
a bijection whose inverse goes as follows. We split each dimer of the pyramid not in the 
right projection into a pair of up- and down-pointing triangles, as shown in Fig. 15- (d). 
Each such triangle has its horizontal edge on a line separating stripes. We also add on the 
right a succession of 6 — a up-pointing triangles adjacent to the lines y = a,a-|-l,...,6— 1. 
For each line y = i, we record the sequence of up- and down-pointing triangles. The walk 
is reconstructed by starting from height a and using successively the triangles in the above 
sequences. More precisely, at each position y = i, the walker makes an up- or down-step 
according to the up-or down-pointing nature of the first unused triangle along the line 
y^i. 

In the above bijection, the descending steps of the walk from height i to height i — 1 
are in one-to-one correspondence with the dimers in the stripe i, except for those of the 
right projection. Note that these descending steps may correspond either to the dropping 
or to the pick-up of a pebble. In particular, weighting each descent of the walk from height 
i to height i — 1 hj Ri (resp. if z < a) amounts to weight each dimer in the stripe i by 
Ri (resp. if z < a), except for those dimers in the right projection of the pyramid (resp. 
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half- pyramid). This immediately yields the relations (A.l). 



Appendix B. Useful identities 

We have the following useful identities, proved by arguments similar to that of Sect. 3. 2. 
For j > 1 and k >0: 

j k 
y^(-iy-'-Un-j,n+j-l{j -i){ Y\ -Rn+fc+1-2^ )-^n-fc,n+fe-l(2i - 1) 

i=0 e=i 

= ^n-j,n+j-l{j) X Sj>k X 5j=k mod 2 

We immediately deduce from this identity the inversion relation 

3 

J2i-'^y~''^n-j,n+j-l{j - i)(Z^-fc,n+fc-l(2z - 1) 

— Rn-k-lRn+k+lZn-k-2,n+k+li'^'i' — 1)) = ^j,k 

The inversion of Eq.(6.1) into Eq.(6.3) follows from this identity. 
For j > 1: 



i=0 



(B.3) 

= —^n-j,n+j-l{j) X Sj=0 mod 2 + JJ^ Rn-£ 



This relation, together with Eq.(B.l) for k = 0, allows to prove Eqs. (6.10) and (6.11) by 

noting that Z„_i,n_i(2z) = Zn,n-l(2z - 1) + -Rn-l-^n-2,n-l(2^ - 1). 
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